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Full resonant transmission of semi-guided planar waves 
through slab waveguide steps at oblique incidence 
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Abstract: Sheets of slab waveguides with sharp corners are investigated. By means of rigorous 
numerical experiments, we look at oblique incidence of semi-guided plane waves. Radiation losses 
vanish beyond a certain critical angle of incidence. One can thus realize lossless propagation through 
90-degree corner configurations, where the remaining guided waves are still subject to pronounced 
reflection and polarization conversion. A system of two corners can be viewed as a structure akin to 
a Fabry-Perot-interferometer. By adjusting the distance between the two partial reflectors, here the 
90-degree corners, one identifies step-like configurations that transmit the semi-guided plane waves 
without radiation losses, and virtually without reflections. Simulations of semi-guided beams with 
in-plane wide Gaussian profiles show that the effect survives in a true 3-D framework. 

Keywords: integrated optics, slab waveguide discontinuities, thin-film transitions, 90-degree wave¬ 
guide corners, vertical couplers, numerical/analytical modeling. 
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1 Introduction 


In a conventional 2-D setting, any abrupt discontinuities in a high-contrast dielectric optical slab waveguide 
typically lead to pronounced radiative losses. This applies also to waveguides with sharp corners. Smooth 
transitions in the form of waveguide bends (B E help to reduce the losses, but also increase the size of the 
structures. Perhaps this is why the intriguing propagation characteristics 01 of defect waveguides in photonic 
crystals, and of 90° bends made of these, attracted so much attention, despite their complexity what concerns 
simulation, design, and experimental realization. With this paper we intend to show that very similar effects can 
be achieved by much more modest means. To this end, corner and step-like structures as illustrated in Figure Q] 
are considered, for out-of-plane guided, in-plane unguided plane waves at oblique angles of incidence. 



Figure 1: Oblique incidence of vertically 
guided, laterally unguided plane waves on 
a 90° corner in a slab (a), and on a step- 
discontinuity consisting of two of these cor¬ 
ners (b); incidence at angle 0. These are 
meant to be Ey-Ez-infinite slabs, folded once 
(a) or twice (b). 


Section E introduces a rather general setting of a slab waveguide “discontinuity”, comprising an interior region 
with in principle arbitrary permittivity, that connects half-infinite slab waveguides at arbitrary positions and 
angles, but aligned such that the entire structure is homogeneous along one lateral coordinate (the y- axis in 
Figure Q]). A general form of Snell’s law applies to pairs of slab modes supported by these access channels. 
Depending on the modal properties of the incoming and outgoing slabs, one identifies critical angles of inci¬ 
dence that border on regimes with vanishing radiative losses, or on regimes with single-polarization or single 
mode propagation. Constituting the basis for early concepts of integrated optics 021, this effect has been 
known for more than four decades (cf. e.g. Refs. e in mu uni), but, to the best of our knowledge, has never 
been applied to configurations other than simple waveguide facets or transitions between co-aligned slabs with 
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different cross sections, in particular not to configurations with orthogonal access slabs, such as for the present 
corners or steps. 

In Sections [3114| we specialize this to high contrast silicon / silicon-oxide slabs. The numerical simulations rely 
on the vectorial implementation of a scheme for 2-D quadridirectional eigenmode propagation (vQUEP) as 
described in Ref. m, based on a former scalar QUEP-variant ma, and building on the bidirectional concepts 
of Refs. mm. To investigate in how far the concepts extend to a more practical 3-D setting with at least 
weak lateral confinement, the propagation of semi-guided, laterally wide Gaussian wave bundles is considered 
in Section \5\ The vQUEP solver GU has been extended accordingly. A preliminary account of the present 
findings, highlighting the mechanism for full-transmission across the step, has been given in Ref. lfT31l . 

For the present paper our primary interest is in the fundamental effect itself. However, structures as discussed 
here could e.g. enable power transfer between layers at different vertical levels in an 3-D integrated optical 
environment, e.g. in a framework of silicon photonics [14, 15]. Beyond the standard evanescent wave interac¬ 
tion ED between the layers, leading to devices measured in centimeters, for vertical distances below 250 nm 
[H6lL or to shorter devices, but for a vertical separation of merely a few tens of nanometers iT7L other concepts 
for vertical coupling include specifically tapered core shapes Ifl8lfl9ll . radiative power transfer through grating 
couplers [20], or even resonant interaction through vertically stacked microrings [2JJ. We believe that the con¬ 
cepts from the present paper could offer here a viable, simple and robust alternative, in particular for situations 
where large lateral extensions of optical channels can be afforded, or if entire layers at different levels need to 
be uniformly “flooded” by light, e.g. for purposes of optical pumping of active devices. 


2 Slab waveguide discontinuities at oblique incidence 

We start with an abstract look at a “discontinuity” in a slab waveguide, as introduced in Figure |2j The structure 
comprises the slab that supports the incoming and possible reflected waves, possible further slabs that support 
outgoing waves, and a central region with — what concerns the arguments in these paragraphs — arbitrary 
properties, hinted at by the dark patch in panel 0a). In case of the configurations of Figure Q] that region covers 
the actual corner, or the corners and the vertical segment, respectively. The entire structure is supposed to be 
constant along the lateral y- axis; we choose the z-direction as the axis of the incoming waveguide. 





Figure 2: Oblique incidence of semi-guided plane waves on a general “discontinuity”, the darker central region; (a) cross- 
section (side-) view, (b) top view of the input waveguide, and (c) generic “top view” of any of the output slabs. Cartesian 
coordinates x, y, z are oriented such that the incoming wave propagates with effective mode index Nj n at an angle 6 in the 
horizontal y-z-plane, with a mode profile that depends on the vertical x-coordinate. Outgoing waves with effective mode 
index N out are observed at angles 6 0U t propagating in the y-(-planes, where ( represents a local coordinate along the slab 
core in the x-z-cross-sectional plane. The structure is constant along the y-axis. 

This concerns time harmonic fields ~ exp(i out) with angular frequency u = kc = 27tc/A, for vacuum 
wavenumber k , vacuum speed of light c, and vacuum wavelength A. For the incoming semi-guided wave 
we assume a polarized guided mode with vectorial profile 'f'[ n (k y ;x) (see Refs. 1171 [TITO and effective mode 
index N- m , supported by the incoming slab, propagating at an angle 6 in the y-z- plane (cf. Figure Ob), also 
Figures [HO 0) . The incoming wave thus relates to a field dependence 

- ^ in (k y ,x) exp(-i [k y y + k z z)), (1) 
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with k y = kN[ n sin A, k z = kN m cos 9 , and k 2 + k 2 = k 2 N^. Due to the homogeneity of the problem along y , 
this harmonic ^-dependence applies to all electromagnetic fields, at all positions; i.e. the global solution can be 
restricted to a single spatial Fourier component with wavenumber k y , here given by the angle of incidence 9. 

Next we look at a particular outgoing mode with profile \Fout and effective index A^ ou t- As hinted at in 
Figure (2a), this can be an actual guided mode supported by one of the slabs with coordinates £i, £ 2 , or 
£4 = —z, but just as well some non-guided, radiated wave propagating in the y-^ 3 - plane. The outgoing field 
(cf. Figure (2c)) is of the form 

~ ^out (ky, .) exp(-i (k y y + &$£))> ( 2 ) 

where the wave equation lITTl requires k 2 N 2 ut = k 2 + k 2 for the cross-sectional wavenumber k^, still with 
k y = kN[ n sin 9 fixed by the incoming field. 

Depending on the angle of incidence, and specifically for each outgoing mode, one thus has to distinguish 
between two cases. If the effective mode index 7V 0Ut is sufficiently large, i.e. if k 2 N 2 ut > k 2 , the outgoing field 
propagates at an angle 9 0 ut with wavenumber = kN 0Ut cos A out , where the effective indices and the angles of 
incidence and refraction are related by Snell’s law: 

A out sin flout = Ain sin A. (3) 

Depending on the properties of the supporting slabs, outgoing waves are thus observed each at its own specific 
angle. 

If, on the other hand, the effective mode index N out is smaller, i.e. if k 2 N 2 ut < k 2 , the outgoing field becomes 
evanescent with an imaginary wavenumber k £ = —i ^Jk* 2 — k 2 N 2 ut . For the cross-sectional problem, these 
waves decay with growing distance £. In particular, these evanescent outgoing fields do not carry optical 
power ll 22 ll (but note that they contribute significantly to the total field around the discontinuity). 

An increase of the angle of incidence A, starting from normal incidence A = 0, can thus cause a change of a 
mode’s type from ^-propagating to ^-evanescent. For an outgoing mode with effective index N out < N- m , this 
happens if k 2 N 2 ui — k 2 N 2 n sin 2 A. Hence, for given input field and for each individual outgoing mode, by 

sin A cr = A out /A in (4) 

one can define a characteristic angle A cr , such that this outgoing mode does not carry power for incidence at 
A > A cr beyond that critical angle. 

With a view to the corners and steps of Figure [H we now consider a structure where simple symmetric three- 
layer dielectric slabs, with core and cladding refractive indices n g and rib, constitute the (identical) access 
waveguides (cf. Figures [2 0). The dimensions are assumed to be such that merely the fundamental modes of 
both polarizations, with effective indices Ateo and Atmo, are guided, where n g > Ateo > Atmo > ^b- The 
structure is excited by the fundamental TE mode. In line with the former arguments, one can then conclude: 

• All modes of these waveguides that relate to radiative waves, with oscillatory behaviour in the cladding 
regions (“cladding modes”), have effective indices below the upper limit rib of the radiation continuum. 
Their characteristic angles (0) are smaller than the critical angle flb, defined by sin Ab = nb/A TE o, asso¬ 
ciated with the background refractive index. Consequently all radiation losses vanish for incidence at 
angles 9 > 9^. 

• All TM polarized modes supported by these waveguides have effective mode indices below Atmo- Their 
characteristic angles © are thus smaller than the critical angle fl m , defined by sinA m = Atmo /Ateo , 
associated with the fundamental TM wave. Consequently, for incidence at angles 9 > 9 m , all incoming 
optical power is carried away by outgoing fundamental guided TE waves. 

Note that these arguments, as well as the former less specific statements, rely solely on the modal properties 
(OD), © of the access waveguides, i.e. the reasoning applies to configurations with — in principle — arbitrary 
interior and arbitrary extension of the region that connects the slab waveguide outlets. 
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2.1 Formal problem 


We complement these more general considerations with a brief look at the rigorous equations, as discussed in 
Refs. mm. The problem is governed by the Maxwell curl equations in the frequency domain for the electric 
field E and magnetic field H , for uncharged dielectric, nonmagnetic linear media with relative permittivity 
e = n 2 , for vacuum permittivity eo and permeability /iq. 

curl E = —iujfioH, curl H = icceeo E. (5) 


The properties of a //-homogeneous structure d y e = 0, and of a corresponding harmonic //-dependence 



(x,y,z) 



(x,z) exp(-i k y y) 


( 6 ) 


of all fields, with the wavenumber k y = kN m sin 9 given by the incident wave, then lead to a system of vectorial 
equations on the x-z-cross section plane, here formulated for the two “transverse” electric field components, 


d x -d z e - d z d x ^ 


d x -d x e + d 2 z 

6 1 

d z -d x e-d x d z d 2 + d z -d z e 


+ k C e ff 


= 0 , 


(7) 


with an effective permittivity e e ff depending on the angle of incidence: 


e e ff(x, z) = e(x, z) - N 2 n sin 2 9. (8) 

Eqs. ©, ® are to be solved on a 2-D domain, with boundary conditions that are transparent for outgoing 
guided and nonguided waves, and that can accommodate the given influx. Note that the problem coincides 
formally with the equations for the modes of channel waveguides with 2-D cross sections. The wavenumber k y 
appears in place of the channel mode propagation constant. Here, however, the equations need to be solved for 
the nonstandard boundary conditions, and need to be treated as a nonhomogeneous problem with the influx as 
a right-hand side, not as an eigenvalue problem as in the case of channel mode analysis. 

Details of the quasi-analytical solver (vQUEP, 112411 ). that has been applied to generate the numerical results of 
this paper, can be found in Refs. lfl2l IT]. In all cases the numerical parameters have been selected such that 
results are converged on the scale of the figures as given, where the power balance (conservation of energy) 
serves as one of the criteria for convergence. With the exception of the configurations with pronounced radiative 
losses (small angles of incidence in Figure |4]), relatively tight computational windows are sufficient, still with 
transparent-influx-boundary-conditions that incorporate the incoming waves (other non-guided, radiated fields 
are suppressed). Similar to guided mode analysis, also here the optical fields are well confined around the 
guiding cores. 

The vanishing of radiative losses for wave incidence beyond the critical angle can be understood just as well in 
terms of these formal equations l23l . In regions with local constant permittivity d x e = d z e — 0, Eqs. © reduce 
to the scalar Helmholtz equation (d 2 + d 2 ) (j) + k 2 e^cj) = 0, valid separately for all components 0 = Ej,Hj of 
the optical electromagnetic fields, with the angle-dependent effective permittivity ® . Wave incidence beyond 
the critical angle 9 b then leads to a negative local effective permittivity, and consequently to evanescent wave 
propagation, in the cladding regions with refractive index n 

The problem © should be distinguished from the standard 2-D Helmholtz problems for TE and TM waves, 
where the global solutions can be represented by a principal scalar field. The latter emerge from Eqs. © if 
both the structure and the solutions are assumed to be constant along one coordinate axis. In our case that 
corresponds to normal wave incidence 9 — 0. Still the familiar scalar TE- and TM-slab modes constitute local 
solutions of Eq. ©, if their vectorial mode profiles are properly rotated 01j], and if their cross sectional 
propagation constants k z , k £ are modified according to the global wavenumber k y . 
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3 Corner discontinuity 


Next we specialize to the corner configurations. Figure |3] introduces a set of parameters that relate to silicon 
thin film cores in a silicon-oxide background at a typical near infrared wavelength. 



Figure 3: Slab waveguide with a straight fold, incidence of ver¬ 
tically (x) guided, laterally (y, z) unguided plane waves at angle 
0, cross-section (a) and top view (b). Parameters: incoming TE- 
polarized light at vacuum wavelength A = 1.55 pm, refractive in¬ 
dices rig = 3.4 (guiding regions), = 1.45 (background), core 
thickness d — 0.25 pm. 


Figures @1 El and Table Q] collect the results of a series of rigorous vQUEP simulations for the corner structures. 
Figure 0] shows the variation of the reflectances Rje, Rjm and transmittances Tte, Ttm of the fundamental 
polarized modes with the angle of incidence. As to be expected, only moderate levels of reflection and trans¬ 
mission (cf. Table |T]) are observed at normal incidence 6 = 0, with most of the optical power being lost to 
radiated fields. Pronounced radiative losses are also present at larger angles below 9\>, where similar amounts 
of power are carried upwards by TE- and TM-polarized waves. 



Figure 4: For the corner of Figure 0 relative reflected (R, dashed) 
and upwards transmitted (T, solid) guided power carried by TE (upper 
panel) and TMpolarized waves (lowerpanel), as a function of the angle 
of incidence 6. The critical angles 6 5 = 30.45° and 6 m = 51.14° are 
indicated. 


As predicted in Section 0 radiation losses vanish entirely Rje + Rtm + Tte + Ttm = 1 for 6 > 9\>, beyond 
the critical angle 9\> for propagating waves in the cladding. Power conversion to outgoing TM waves reaches a 
strong maximum in between 9\> and 9 m , and is suppressed entirely beyond the critical angle 9 m for propagating 
TM waves, where all outgoing waves are purely TE polarized: Rje + Tte = 1 for 9 > 9 m . For even higher 
angles of incidence one observes a maximum in the TE transmission, which vanishes for grazing incidence as 
9 approaches 90°. 



Figure 5: Field profiles \E\ for the corner discontinuity of Figured for angles of incidence 0 = 0°, 26°, 41°, 68 °. 

Figure El provides a few examples of fields observed for these typical cases. The plots show absolute values 
of the electric field vector \E\, with the color scales adjusted such that variations in regions with small field 
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strengths become visible. Note that, irrespective of the critical angles, these rigorous solutions contain con¬ 
tributions from the complete sets of local vectorial slab modes, of both polarizations, which are used as the 
expansions bases |[TT1 . If beyond a related critical angle, however, the then evanescent modes stay localized 
around the corner region, not contributing to the “far-fields” in the outlets. 

While pronounced radiative losses are evident in parts (a) and (b) for 6 < Q\>, all fields remain confined around 
the core regions in (c) and (d) for angles 6 > Q\>. Outgoing TM waves lead to strong electric fields immediately 
outside the cores of the vertical channel (“highlighted” core edges) in panels (b) and (c), but are absent for the 
scalar problem at normal incidence (a) and for 6 > 6 m in (d). 

Tabled] compares values for reflectances and transmittances, together with selected parameters, for the config¬ 
urations of all field plots in this paper. Although here we consider the corners merely as building blocks for the 
steps of the next sections, a lossless corner structure as the one of Figure ®c), that channels a total of 74% of 
the horizontal input into the vertical, might be of practical interest in its own. 
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Table 1: Reflectance and transmittance values Rje, Rtm , Tte> Tjm for the polarized fundamental slab modes/ for wave 
bundles made of these, and varying parameters (angle of incidence 0, vertical layer distance h, bundle widths W y , W cr 
along y and along the beam cross section) for the configurations relating to the field plots in this paper. 


4 Step discontinuity 

Combination of two of the former corners gives a step structure that connects optical layers at different ele¬ 
vations. Figure [6] introduces coordinates and parameters. Aiming at efficient power transfer, we focus on two 
of the lossless cases from Section [3} for 0 = 41° (Figure 0c), this will be named configuration CO) and for 
6 = 68° (Figure0d), labeled configuration Cl). Both relate to transmission maxima in Figure® For each case 
we assume that the two combined corners are identical. 


Figure 6: A step configuration, cross-section (a) and 
top view (b). Parameters are as given for Figured 
with an additional value h for the vertical core dis¬ 
tance. 

With the vertical distance h between the horizontal layers, only one new parameter appears. Parts (a) and (c) of 
Figure [7] show the transmittance and reflectance of the step as a function of this distance, calculated by rigorous 
vQUEP simulations for the full step structure. For both configurations, the angles of incidence exceed the 
critical angle 6 \>, hence the steps are lossless. 

We look at configuration (Cl) first. Since the related corner structures (cf. Figure ®d), entry in Tabled]) trans¬ 
mit and reflect only TE waves, only the upward and downward propagating versions of the fundamental TE 
mode mediate between the two corners in the vertical part of the step. One might thus compare the step with 
a standard Fabry-Perot interferometer [25 ], where the corners play the role of the (identical) partial reflectors. 
Consequently, the scan over the vertical separation h reveals multiple-beam fringes resulting from the interfer¬ 
ence of the repeatedly up- and downwards reflected TE modes. Maxima with virtually full transmission appear 
for certain equidistant heights. 
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Figure 7: For the steps of Figured TE-input: polarized modal transmittance (Tje, Tjm) and reflectance (Rje, Rtm) 
versus the step height h (a, c) for given incidence angle 6, and versus the angle of incidence for given step height (b, d). 
Thin dash-dotted lines indicate the values 6 = 41°, 68° and h = 1.83 /iffl, 2.15 pm selected for the configurations (CO) 
and (Cl) of Figure^ 


Figure [8fb) shows the field pattern for an optimum configuration with full TE-to-TE transmission over a vertical 
distance of h = 2.15 /im. The field is confined around the core regions, without any beating pattern, i.e. with 
purely forward waves, in the horizontal segments, and with strong resonant, mostly standing waves in the 
vertical slab. 



- 3 - 2-1 0 1 2 3 - 3 - 2-1 0 1 2 3 

(a) z [|im] (b) z [jum] 


Figure 8: Field profiles \E\ for two step configurations as in Figure\6 1 with (almost) full transmission. For the fundamental 
guided slab modes one observes transmittance and reflectance levels as given in Table{l\ 

Regarding configuration (CO), at 0 — 41° and for TE input, the corners transmit the fundamental modes of 
both polarizations (with a dominant TM part, cf. the entry of Figure [5^c) in Table Q]). Hence, in a step structure 
made of these corners, one must expect that TM as well as TE waves play a role in the interference in the 
vertical segment. Accordingly, the scan over the vertical separation in Figure [8^a) shows a much less regular 
dependence. Still, for a distance h = 1.83 /xm, a lossless configuration with virtually full transmission can be 
identified. The transmission is mostly TE polarized (reciprocity arguments apply LIT]), with a TM contribution 
of about 2%. The field plot in Figure [8ja) shows a major TM contribution to the resonance in the vertical slab, 
visible through the strong electric field contributions next to the vertical core edges, similar to Figure [5fc). 

In principle, given the properties of the constituting corners in the form of scattering matrices, optical transmis¬ 
sion through the step structures can be analyzed analytically. For configuration (Cl), where propagating TM 
waves are suppressed, the corners are represented by a 2 x 2 matrix that relates to incoming- and outgoing TE 
waves. Configuration (CO) requires a 4 x 4 matrix to accommodate TE as well as TM waves. The model would 
be applicable to steps of sufficient height, where evanescent fields around the corner regions do not play a role. 

As a resonant effect, the property of full transmission must be expected to be sensitive, to some degree, to all 
parameters that enter. With a view to Section[5]we select the input angle for a further parameter scan. According 
to panels (b) and (d) of Figure [71 although less perfect in peak performance, configuration (CO) might turn out 
to be more robust than (Cl), here concerning variations of the angle of incidence. 
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5 Semi-guided beams 


All examples discussed in Sections [3] and 0 concern semi-guided plane waves that extend infinitely in the ±y 
direction. We now look at bundles of the former solutions, for a typically small range of angles of incidence, 
or of wavenumbers k y , respectively, with the aim of modeling the incidence of semi-guided, laterally wide but 
localized beams on the corners and steps. These are wave packets of the form 


E 

H 


J (x, y,z) = A j exp 


{ky kyo)* 


wt 


• ^ in (k y ;x)exp(-ik z (ky)(z-z 0 )) +p(k y ]x,z)^ exp(-ik y (y-y 0 )) d k y , (9) 

with a Gaussian weighting of half width centered around a primary wavenumber k y o = kN[ n sin 9 q, for 
primary angle of incidence 9q, and an (arbitrary) amplitude A. Phase factors have been introduced to position 
the focus in the y-z -plane at (yo, zo). The central term in curly brackets represents the numerical 2-D vQUEP 
solution, formally separated into a contribution of the incoming field with mode profile \l/i n , and a remainder 
p. The integrals are evaluated by numerical quadrature 126ft . We thus arrive at approximations of ^-localized 
true 3-D solutions, as a basis for the results in Figures l9l and flOl 



Figure 9: Optical electromagnetic energy density for incidence of Gaussian wave bundles (Eq. (0) with a cross section 
width W cr = 10 fim on a corner discontinuity as in Figure\3\ for principal angles of incidence (a) 9q = 41° and (b) 
9 q = 68°. The panels show cross section views (i), bottom views (ii), and views of the plane through the center of the 
vertical slab (Hi). Further parameters, and transmittance and reflectance levels for the bundles are given in Table\7} 


For a proper choice of wave packet parameters it is instrumental to evaluate the incident part (everything except 
p) of Eq. (|9]) a little further. Assuming a small spectral width Wk, and neglecting the effect of the rotation of 
the vectorial mode profile (by replacing \Pi n (fc y ; x) by i &m(kyO m i x ))> the incident field can be expressed as 


8 












Figure 10: Optical electromagnetic energy density for incidence of Gaussian wave bundles © with cross section widths 
ofW cr = 10 gm (COn, a), W cr = 45 gm (COw, b), W cr = 10 gm (Cln, c), W cr = 180 gm (Clw, d), on step discontinuities 
as in Figure^ for principal angles of incidence 0 q = 41° (CO, a, b) and 0 q = 68° (Cl, c, d). The panels show cross 
section views (i), views of the plane through the center of the vertical segment (ii), and bottom views at the levels of the 
lower (Hi) and of the upper slab (iv). Modal transmittances and reflectances for the bundles are given in Table\l J 
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( 10 ) 



^2 A 

*Wy' a * 




{y ~ yo) - kf 0 ( z - z °)) ^ 

(Wy/ 2) 2 


• ^inOyo; z) exp (-i(feyo(y - yo) + Z^oO - zo))), 


with k z o = kN[ n cos 0 q* In the y-z- plane this is a Gaussian beam with fields of full width W y along y at 1/e- 
level, where Wy = 4/w&. By introducing the cross section position c and longitudinal position l , relative to the 
focus, as new coordinates through y = yo + Z sin + ccos ^ + l cos #o — csin one can write the 

incoming field in the more succinct form 

(h). (^>^0 ~ exp ^~ (-yy C ^ 2 ) 2 ^ &in( k yo;x) exp (~ikN in l). (11) 

Here the cross-section-width lW cr (full width of the field at 1/e-level, at focus, in the direction perpendicular 
to the beam axis) is related to the y- width by W cr = W y cos 6% Both quantities might be practically relevant, 
hence both values are listed in Tabled] 

Figures l9l and ITOl collect corresponding results for our corner and step structures. The plots show the electro¬ 
magnetic energy density (eoe\E\ 2 + yo\H\ 2 )/4: as a physically more relevant quantity. Uniform color scales 
have been adopted for all three or four panels that relate to different views of the same configuration. In some 
cases the strong resonant fields visible in the longitudinal views lead to rather dark incoming/outgoing beams 
(although these carry power of unit order). The single contour at 2% of the energy density maximum is meant 
to better accentuate the form of these beams. 

According to the values in Table \T\ both corner configurations work almost identically for bundles with cross 
section width W cr = 10 ym and for incoming plane waves. The reflected and transmitted beams in Figure [9] 
remain nicely confined. This can be explained by the mere weak angular dependence of the transmission curves 
in Figure 0] around the maxima selected for the corner configurations (CO) and (Cl). 

In contrast, for the same wave bundles of width lW cr = 10 ym, the performance of the step structures de¬ 
teriorates, when compared to the case of plane wave incidence. Here the strong angular dependence of the 
resonant transmission maxima in Figure [7] becomes relevant. The reflected beams in parts (a, iii) and (c, iii) of 
Figure [TO] show pronounced sidelobes. Merely the central part of the incoming bundle appears to be transmit¬ 
ted. The effect is more pronounced for configuration (Cl) than for (CO) due to the narrower transmission peak 
in Figure |7{d). 

Hence we choose wider input beams for both configurations as our last examples. The cross section widths 
of Wcr = 45 ym (CO) and W cr = 180 ym selected for panels IToTb, d) roughly correspond to an angular 
range around 9o with TE-transmission above 90%, i.e. and W y , W cr , respectively, are adjusted such that 
T te (9) > 90% for all 9 <G [0o — 9 q + A0], for A 9 = w^/{kN m cos 0o). According to Tabled] for incoming 
beams this wide, the step configurations come reasonably close to their ideal plane-wave performance. 


6 Concluding remarks 

Semi-guided plane waves at sufficiently high angles of incidence propagate across arbitrary straight slab wave¬ 
guide discontinuities without any radiation losses. This effect can be understood in terms of a variant of Snell’s 
law, that relates effective mode indices and angles of incidence and refraction in pairs of access channels. For a 
quantitative analysis of these configurations, the frequency-domain Maxwell equations reduce to a vectorial 2- 
D system, to be solved on a computational window with transparent-influx boundary conditions. Our rigorous 
quasi-analytical solver enables the convenient analysis of structures with general rectangular permittivity distri¬ 
butions. Simulations of high-contrast Si/Si02 corners lead to the identification of two example configurations 
with local transmission maxima. A Fabry-Perot-like resonance effect permits to combine two identical corners 
into step configurations with full transmission of the incoming semi-guided plane wave, i.e. with (numerically) 
ideal behaviour. Bundles of the former solutions can serve as examples of what happens to semi-guided, lat¬ 
erally localized beams, when incident on the corners or step structures. Our examples show that transmission 
properties very close to the laterally infinite case can be realized with beams of sufficient width. 
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Once a working step configuration has been identified, extension to further, perhaps also intriguing examples 
is obvious. Due to the waveguide symmetry, the transmission properties of the corners is not affected if one 
changes their direction. This enables u-tum-like structures as in Figure \TT\a). Resonant vertical segments in 
the form of steps or u-turns can be connected by pieces of horizontal slabs of arbitrary length. One is thus led 
to the bridge-, staircase-, or s-bend-like configurations of Figure ITTTb. c, d). Our vQUEP simulations predict 
transmittance levels T > 99% in all cases. 



Figure 11: Field profiles 
\E\ for u-turn, bridge, 
stair, and s-bend -like 
structures, each made of 
double-corners with pa¬ 
rameters as for configura¬ 
tion (Cl) of Section^ 

(0 = 68°, h = 2.15 pm), 
with an (arbitrary) z- 
distance of 3 fim between 
the vertical sections in (b, 
c, d). 


With the exception of the height scan of Section 01 design of the step structures is straightforward, without the 
necessity of further tuning of parameters. Still there is ample room for optimization: Apart from variations of 
layer thicknesses, other shapes could be investigated for the refractive index profile of the corner regions. One 
might aim at a higher transmittance of the single corners, and consequently at a wider angular range of high 
transmittance for the steps with higher transmission of narrower bundles, preferably for single polarization 
configurations such as the present (Cl). One might also consider structures with true lateral guiding. This 
would require the preparation of wide channels with weak lateral contrast (e.g. by shallow etching, indiffusion, 
direct laser writing, or other processes) along the path of the present beams. 
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